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Abstract 

In the Minimal Supersymmetric Standard Model (MSSM), the bottom Yukawa coupling of the 
Higgs boson can considerably deviate from its Standard Model prediction due to non-decoupling 
effects. We point out that the ratio of the Higgs boson decay branching fraction to a bottom quark 
pair and that to a lY-boson pair from the same production channel is particularly sensitive to 
large additional MSSM Higgs boson mass regions at future electron-positron colliders. Based on 
this precision measurement, we explicitly show the indirect discovery reach of the additional Higgs 
bosons according to planned programs of the International Linear Collider. 
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In July 2012, discovery of a new particle with a mass of 125 GeV was announced by the 


ATLAS and CMS collaborations at the CERN Large Hadron Collider (LHC) [l|. Properties 
of the discovered particle have turned out to be consistent with those of the Higgs boson of 
the Standard Model (SM) within the error. The SM has been experimentally conhrmed as 
a low-energy effective theory that consistently describes every phenomenon at energy scales 
below 0(100) GeV. Although the SM is very successful, several phenomena which cannot 
be explained in the framework of the SM have been known, such as neutrino oscillations, 
the abundance of dark matter and baryon asymmetry of the Universe. In order to account 
for such strong evidence we need to go beyond the SM. 

It should be emphasized that many new physics models which can explain the new phe¬ 
nomena mentioned above demand extension of the Higgs sector. Given such a non-minimal 
Higgs sector, a SM-like Higgs boson with a mass of 125 GeV is not exactly the SM one, and 
has different characteristic properties from their SM predictions. Therefore, it is plausible 
that hints of new physics are obtained by investigating properties of the discovered Higgs 
boson. In this viewpoint, the discovered Higgs boson is a window to new physics. 

Indeed, one of the best strategies for exploring new physics is to investigate its effects that 
appear in the couplings of the discovered Higgs boson to SM particles. Since the magnitudes 
of these Higgs boson couplings depend strongly on masses and couplings of particles yet to 
be discovered, hngerprinting of the Higgs boson couplings is useful in determining the energy 
scale of new physics and further in distinguishing new physics models behind. At the LHG, 
deviations of the Higgs boson couplings from their SM predictions are currently constrained 
typically with 10% |2|. At future electron-positron colliders such as the International Linear 
Gollider (ILG) j^, the Gompact Linear Gollider (GLIG) jdl and the Future Gircular Gollier 


of electrons and positrons (FGG-ee) , accuracies for measurements of various Higgs boson 
couplings will be signihcantly improved to a percent level or better jo-S]. Loop corrections 
from new particles also give characteristic deviation patterns in these couplings. Preci¬ 
sion measurements at future collider experiments can differentiate such small deviations. 
Therefore, it is crucial to perform hngerprinting of the deviation pattern of the Higgs boson 
couplings based on precise computations including radiative corrections in beyond-the-SM 
(BSM) models for estimating the scale of new physics indirectly. 

Among ever proposed BSM models, supersymmetric (SUSY) extension of the SM is one 
of the excellent paradigms j9|. Introducing superparticles whose spins are different from the 
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SM counterparts by one half, an unnatural cancellation between the bare mass squared of the 
Higgs boson and its quadratic ultraviolet divergences is avoided. By imposing the i?-parity, 
the lightest SUSY particle (LSP) is stabilized and thus can be an excellent candidate for dark 
matter. The SU{3)c, SU{2)i and U{1)y gauge coupling constants can be naturally unihed 
at a high energy scale, inspiring many ideas about supersymmetric grand unihed theories. 
One of the intriguing features is that the Higgs sector of SUSY models must be extended to 
cancel the gauge anomaly. Even in the Minimal Supersymmetric Standard Model (MSSM), 
two Higgs doublets H 2 with opposite hypercharges are required. As a consequence, the 
coupling strengths of the discovered Higgs boson with a mass of 125 GeV to SM particles 
deviate from their SM predictions, and have the same structure as in the type-H two Higgs 
doublet model at the tree level. At the loop level, however, superparticles radiatively affect 
the Higgs boson coupling strengths, giving rise to characteristic deviations patterns different 
from other models. In particular, it has been known that signihcant deviations of the bottom 
and tau Yukawa couplings induced by non-holomorphic radiative corrections do not decouple 
even in the large superparticle mass limit lOHlOj. Due to the non-decoupling effects, the 
indirect reach of the additional Higgs boson mass scale through precision measurements of 
the bottom/tau Yukawa coupling can be much higher than the direct reach at the LHC and 
its luminosity-upgraded version. 

In this letter, we reanalyze the deviation of the Higgs boson couplings enhanced by 


the loop-induced non-decoup 


recent relevant works 


18 


ing effects from a difference perspective. In sharp contrast to 


19|, we utilize the ratio of Br(h —)■ hh) to Br(h WW) from 


the same production channel e“*'e“ vvh as a key observable. Since this combination 
cancels out the uncertainties related to the Higgs boson production cross section and total 
decay rate, the indirect additional Higgs boson reach obtained through this method is more 
expanded than using single coupling accuracies. Since the cross section times branching 
ratios of these decay modes are particularly accurately measurable, our quick method is not 
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inferior to global analysis. As an example, we assume expected accuracies at the ILC 
whose running programs are under discussion. Our goal is to explicitly show the indirect 
discovery reach of the additional MSSM Higgs bosons through the measurement of the ratio 
Br(h —)■ 56)/Br(h WW) at several planned ILC programs. Such an analysis plays an 
important role in properly assessing the capability of future electron-positron colliders for 
exploring new physics effects. 
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We start briefly reviewing the MSSM Higgs sector, which consists of two Higgs doublet 
scalars. The electroweak symmetry is broken when the neutral components of and H 2 
develops vacuum expectation values (VEVs). Down-type quarks and charged leptons (up- 
type quarks) acquire their masses from the VEV (H^) = Vi/\/2 ~ The ratio 

tan/3 = V 2 lv\ is an important quantity on which many observables depend, and the sum 
of the squares is given by = v\ ^ v\ = (246 GeV)^. After the electroweak symmetry 
breaking, there are two GP-even Higgs bosons h and H, a GP-odd Higgs boson A and 
charged Higgs bosons as physical particles. The lighter GP-even Higgs boson h is 
identified as the discovered SM-like Higgs boson at the LHC. At the tree level, its predicted 
mass is smaller than the mass of the Z-boson as the Higgs self-interactions are given by 
SU{2)l and P(l)v gauge couplings. However, the radiative correction signihcantly lifts 
the mass of the Higgs boson and can account for its observed mass of 125 GeV although 
superparticles are required to have a few TeV masses 20|]. The mixing angle a between 
the two GP-even Higgs bosons is also affected by the radiative correction. For detailed 
discussion about the radiative corrections to the masses of the GP-even Higgs bosons and 
its mixing angle, see, for example. Ref. 2l|. By virtue of supersymmetry, H, A and are 
degenerate in mass in the large mass limit. We consider the GP-odd Higgs boson mass 
as a measure of the mass scale of these heavy MSSM Higgs bosons. 

Let us discuss non-decoupling effects of radiative corrections to Yukawa couplings in the 
MSSM. In this letter we focus on the bottom Yukawa coupling because it receives large 
radiative corrections from the strong coupling and the top Yukawa coupling as well as 
future linear colliders are capable of measuring the h ^ bb decay rate with considerable 
accuracy. In SUSY models, the holomorphy restricts possible form of the superpotential. 
As a consequence, the bottom quark does not couple to H 2 at the tree level. However, effects 
of SUSY breaking induce a non-holomorphic bottom Yukawa coupling A A;, radiatively in 
addition to holomorphic radiative corrections 6Xb jl^ : 


— Cb — (Aft -|- 6Xb)bjiHiQ]b + AXbbjiQLH^ , 


( 1 ) 


where b^ and Ql denote the right-handed bottom quark and the third generation quark 
doublet, respectively. Notice that such non-holomorphic couplings are absent in the type-H 
two-Higgs-doublet model due to the hypothetical Z 2 symmetry. The mass of the bottom 
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quark is also shifted by the radiative corrections as 



( 2 ) 


where 



( 3 ) 


The coupling of the SM-like Higgs boson to the bottom quark is obtained as 



— Aft J (1 + cot a cot /S) 


( 4 ) 


Dominant contributions to the non-holomorphic coupling stem from sbottom-gluino and 
stop-chargino loop diagrams and are given by 



where /i, M 3 and At are the higgsino mass, the gluino mass and the trilinear stop A- 
parameter, respectively. The typical mass scale of the superparticles in the loop diagrams is 
collectively denoted by m-gusv Notice that there are two mass scales characterizing effects 
of new physics in the MSSM. One is the mass scale of superparticles msusY, and the other 
is the mass of the OP-odd Higgs boson rriA. In the limit where superparticles are heavy 
fi ~ M 3 ~ ~ msusY S> mz with fixed niA, the bottom Yukawa coupling is approximately 

given by 



( 6 ) 


Therefore, SUSY loop corrections do not decouple for small rriA, and are enhanced for large 
tan 13. On the other hand, in the limit where the mass of the OP-odd Higgs boson is large, 
we obtain 
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and thus SUSY radiative corrections decouple. It should be noticed that two-Higgs-doublet 
models can also produce non-decoupling effect on the Higgs boson couplings, which originates 
from extra Higgs boson loop diagrams 3-26|. 

We here describe the method for our numerical analysis. We utilize the public code Feyn- 
Higgs2.10.3 to compute mass and mixing parameters as well as decay branching fractions 
in the MSSM Higgs sector 27| The masses of the UP-even Higgs bosons and the mixing 
angle are computed at the two-loop level in the on-shell scheme. Branching ratios of the SM- 
like Higgs boson are computed at the one-loop level including QED and QCD corrections. 
We perform a random scan of MSSM parameters avoiding tachyonic scalar masses. To cir¬ 
cumvent large UP-violating processes, we constrain our analysis to real MSSM parameters. 
For simplicity, we assume that soft SUSY breaking H-parameters are given by the product 
of the corresponding Yukawa coupling constant A/ and a common Aq value as Af = XfAo. 
We consider the mass range of electroweak (colored) superparticles up to around 1000 GeV 
(4000 GeV) to cover less tuned parameter sets which evade superparticle searches at the 
future LHG experiment. Since the non-decoupling effects depend on ratios of superparticle 
masses, our results about the Higgs boson couplings are not substantially changed even for 
heavier superparticles. If \At\ is significantly larger than the stop masses, there appear a 
deeper charge/color breaking vacuum than our electroweak vacuum. To avoid complexity of 
the discussion of the meta-stability bound, we constrain |Ho| to be smaller than three times 
the geometric average of the soft stop masses. The scan bounds of the MSSM parameters are 
listed in Tab. [Tl Since radiative corrections from the first two generation sleptons and squarks 
to the Higgs sector are small, we set the soft masses of these sleptons and squarks to 1000 
GeV and 4000 GeV, respectively, instead of performing a random scan. This degeneracy also 
suppresses unwanted flavor chaining neutral currents and lepton flavor violating processes 
adequately. Although the measured mass of the Higgs boson is now fixed to rrih = 125 GeV 
with the error less than 1 GeV 28(], we consider 122 GeV < < 128 GeV as an allowed 

region for its mass. This is a judicious choice to take uncertainties from renormalization 
scheme dependence and higher order contributions to the mass of the Higgs boson into 
account 30|]. From the cosmological viewpoint, it is assumed that the lightest neutralino is 


the LSP and a dark matter candidate. However, the relic density constraint is relaxed in 


^ Changes in the recently updated version FeynHiggs2.11.0 are not relevant to our analysis. 
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TABLE I: Scan bounds on MSSM parameters. 


Parameter 

Scan bounds 

rriA 

[200 GeV, 3000 GeV] 

tan f3 

[1, 60] 


[100 GeV, 1000 GeV] 

"^Q.3,f/3,D3 

[500 GeV, 4000 GeV] 


[0 GeV, 

l/^l 

[100 GeV, 1000 GeV] 

|MJ 

[100 GeV, 1000 GeV] 

1^41 

[100 GeV, 1000 GeV] 

M3 

[1400 GeV, 4000 GeV] 


TABLE 11: Constraints on superparticle masses. 


Observables 

Gonstraints 


[122 GeV, 128 GeV] 

m~ 

> 1400 GeV 

m~o 

X? 

> 90 GeV 

m_± 

Xi 

> 400 GeV 

rrir 

u 

> 640 GeV 

rur 

> 620 GeV 

m~ 

> 90 GeV 


order to allow for any non-standard cosmological scenario. As for the other experimental 


constraints, we impose superparticle mass bounds obtained at the LEP 


31 | and the LHC 


32| . The mass bounds we employ are summarized in Tab. [TTl The flavor violating decays 
Bg —)■ and b ^ s'j mediated by the additional Higgs bosons could also give constraints 

on the SUSY parameter space, as studied in |^. However, such processes can be canceled 
by contributions from squark sector parameters without changing our results. Therefore, 
we neglect these flavor constrains in our analysis. The deviation of the anomalous magnetic 
moment of the muon from its SM prediction can be explained by adjusting the masses of the 
second generation sleptons, which is beyond the scope of this letter. Since our analysis is 
conhned to low-energy SUSY scenarios, we require Xb to be perturbative at the electroweak 
scale. 

Our goal is to show to what extent large scenarios with moderate SUSY parame- 
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TABLE III: Expected accuracy for R at planned stages of the ILC 


CM energy and luminosity 

95%GL 

5a 

= 500 GeV, L = 500 

0.050 

0.125 

Vs = 500 GeV, L = 1600 fb"^ 

0.027 

0.068 

Vs = 1 TeV, L = 2500 fb"^ 

0.021 

0.052 


ter choices can be surveyed through the precision measurements of the cross section times 
branching ratio measurements of the h ^ bb and h —)■ WW modes which follow the pro¬ 
duction vVh. To this end, we introduce the following double ratio of the branching 

fractions oi h ^bb and that of h —?■ WW modes: 


R = 


Br(h ^ WW) 


Br(h —)■ bb) 


Br(h ^ WW) 


MSSM . 


Br(h —)■ bb) 


( 8 ) 


SM 


Its tree level expression is given by 


ll-tati/3cot(;3-a)|2’ 


(9) 


where cot(/3 — a) is given as a function of ttia and mz- As is evident, any deviation from 
the SM prediction R = 1 exhibits the existence of some new physics beyond the SM. By 
taking the ratio of the cross section times branching ratios for bb and WW modes in the 
same production mechanism, uncertainties that originate from the production cross section 
and the total decay rate of the Higgs boson are canceled out. Consequently, constraining 
the model parameter space using the accuracy of the ratio of the decay branching fraction 
R is more powerful than using the converted accuracy of the hbb coupling. The expected 
accuracies for cross section times branching ratio measurements of the SM-like Higgs boson 
through the uuh mode at several ILC stages with the polarization options {Pe-,Pe+) = 
(—0.8,-1-0.3) and {Pe-,Pe+) = (—0.8,-1-0.2) are listed in Tab. 5.4 and 5.5 of Ref. js]. 
Sensitivities for R are computed based on these tables and summarized in Tab. IIIII For 
example, for y/s = 500 GeV and L = 500 fb“^, the accuracy of the ratio at 95% CL is 
estimated as -^/((L7%)M^~(2H%)^ x 2 = 5%, while the accuracy of Br(h —)■ bb) from the 
coupling measurement itself is at best ~ 6.4% ignoring systematic uncertainties (see Table 
6.1 in Ref. [8]). 
















In Fig. [H the indirect 5a discovery reach of the MSSM heavy Higgs bosons through the 
measurement of R at the ILC with the center-of-mass (CM) energy ^/s = 500 GeV and 
the luminosity L = 500 fb~^ is shown in the (m^,tan/9) plane. MSSM scenarios which 
predict more than 5a deviation of R are plotted with the red points. The blue line is the 
corresponding tree level prediction based on i?tree- The indirect reach of the MSSM heavy 
Higgs boson mass can be significantly extended due to the non-decoupling property of the 
non-holomorphic radiative corrections, compared to the tree level results. When the relative 
sign of yU and M 3 (A^) is minus, the deviation of R from unity is strongly enhanced and 
the correction is proportional to tau/S (see Eqs.(j5]) and ([ 6 ])). The solid (black) line shows 
the 95% exclusion limit obtained at the LHC Run-I through direct H/A ^ searches 
. At the LHC Run-H with the luminosity of L = 10 fb“^ accumulated, it is expected 


to exclude up to = 1000 GeV for tan/? = 10 in view of the increase of the production 
cross section 3^. The existence of the heavy Higgs boson can be searched also through 


the H/A —)■ jjAn decay modes at the LHC |35|. Its expected 5a discovery reach with 


L = 3000 fb~^ almost overlaps with the 95% excluded region explored by the r channel. 

The 5a heavy Higgs mass reach expected at the stage with = 500 GeV and L = 1600 
fb“^ and at the stage with ^/s = 1 TeV and L = 2500 fb“^ are shown in Fig. [2] and Fig. |3l 
respectively. These figures show that luminosity-upgraded ILC programs are very powerful 
in measuring the Higgs boson branching fractions and thus sensitive to higher regions 


Before closing, we comment on the parametric uncertainties from the SM parameters 


Theoretical calculation on Br(h —)■ bb) suffers from uncertainties on and a^, and the 
uncertainties on the calculation limit the reach of MSSM parameters by experimental mea¬ 


surements 


37l | . Although current errors on these parameters 


lead to larger uncertainties 


than the expected uncertainties at the future ILC measurements, it is anticipated that fu¬ 
ture developments on higher order calculation and lattice calculation overcome this situation 
until the experimental programs take place 391. 

In this letter, we have presented to what extent the heavy additional MSSM Higgs bosons 
can be indirectly accessible by using precision measurements of the cross section times 
branching ratios for the discovered Higgs boson at the future ILC experiments. We have 
employed the ratio of the Higgs boson decay branching fractions Br(h — bb) /Br(/i —)■ WW) 
as a simple but useful observable for indirectly probing the existence of the heavy Higgs 
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FIG. 1: Indirect 5a discovery reach of the heavy Higgs bosons through the measurement of R at the 
ILC stage with ^/s = 500 GeV and L = 500 fb“^ in the (mA,ta-n/3) plane (red points). For comparison, 
the indirect 5a reach based on the tree-level results is delineated with blue line. The 95% exclusion limit 
obtained at the LHC Run-I through H/A ^ searches is delineated with black line. 

bosons because this ratio is obtained by using the cross section times branching fractions, 
which cancel out the uncertainties of the production cross section and total decay width 
of the Higgs boson. We have computed the ratio R, which can substantially deviate from 
unity even for large rriA regions due to large non-decoupling effects of superparticles, and 
have exhibited the indirect 5cr discovery potential of the additional MSSM Higgs bosons at 
several planned stages of the future ILC experiments. Our analysis shows that the ILC, in 
particular with high luminosity, cover many large rriA scenarios in the MSSM. 

Last but not least, hngerprinting of Higgs boson couplings at the loop level including our 
analysis improves capability of exploring and distinguishing the MSSM from other models 
of the electroweak symmetry breaking, and should be encouraged. 


ILC: Vs = 500 GeV, L = 500 fb '' 
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